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We investigate the steady-state spin diffusion for ultracold spin-orbit coupled 40 1< gas by the 
kinetic spin Bloch equation approach both analytically and numerically. Four configurations, i.e., 
the spin diffusions along two specific directions with the spin polarization perpendicular (transverse 
configuration) and parallel (longitudinal configuration) to the effective Zeeman field are studied. 
It is found that the behaviors of the steady-state spin diffusion for the four configurations are 
very different, which are determined by three characteristic lengths: the mean free path l T , the 
Zeeman oscillation length In and the spin-orbit coupling oscillation length l a . It is analytically 
revealed and numerically confirmed that by tuning the scattering strength, the system can be divided 
into five regimes: I, weak scattering regime ( l T > ln,l a )', II, Zeeman field-dominated moderate 
scattering regime (In <C l T l a ); III, spin-orbit coupling-dominated moderate scattering regime 
(l a <IC l T <S In)’, IV, relatively strong scattering regime (l( : < I T < In, la)’, V, strong scattering 
regime (l T -C In, la, It), with l / representing the crossover length between the relatively strong and 
strong scattering regimes. In different regimes, the behaviors of the spacial evolution of the steady- 
state spin polarization are very rich, showing different dependencies on the scattering strength, 
Zeeman field and spin-orbit coupling strength. The rich behaviors of the spin diffusions in different 
regimes are hard to be understood in the framework of the simple drift-diffusion model or the direct 
inhomogeneous broadening picture in the literature. However, almost all these rich behaviors can be 
well understood from our modified drift-diffusion model and/or modified inhomogeneous broadening 
picture. Specifically, several anomalous features of the spin diffusion are revealed, which are in 
contrast to those obtained from both the simple drift-diffusion model and the direct inhomogeneous 
broadening picture. 

PACS numbers: 67.85.-d, 51.10.+y, 03.75.Ss, 05.30.Fk 


I. INTRODUCTION 

In recent years, spin dynamics including spin relax¬ 
ation and spin diffusion/transport is extensively studied 
in both Boseid— and Fermi^r— cold atoms. For the Bose 
system, the spin dynamics of the Bose-Einstein condensa¬ 
tion has attracted much attention^ 8 For the Fermi cold 
atoms, the systems without^— and with2£3 8 2 spin-orbit 
coupling (SOC) are extensively investigated. In the ab¬ 
sence of the SOC, many interesting phenomena, such as 
the Leggett-Rice effect in unitary gas^^r— and anomalous 
spin segregation in extremely weak scattering limit 
have enriched the understanding of the spin dynamics of 
Fermions. With the synthetic SOC experimentally re¬ 
alized by laser control technique in cold atoms 
the spin relaxation for the Fermi cold atoms with 
SOC has been studied both experimentally^— and 
theoretically—This is partly motivated by the well- 
controlled laser technique, which provides more freedom 
for the cold atoms than the conventional solids. On one 
hand, rich regimes can be realized by tuning the SOC 
strength; on the other hand, not only the interatom in¬ 
teraction can be tuned by the Feshbacli resonance^ but 
also the atom-disorder interaction can be introduced and 
controlled by the speckle laser technique— 

The experimentally realized effective Zeeman field 
and SOC provide an effective magnetic field, which 


reads^ i 26 ' 27 

f2(k) = (Cl, 0, 5 + ak x ). (1) 

In above equation, k = (k x ,k y ,k z ) denotes the center- 
of-mass momentum of the atom; Cl acts as an effec¬ 
tive Zeeman field along the x-direction; S is the Ra¬ 
man detuning, which is set to be zero in our work; 
fl z (k) = ak x represents the k-dependent effective mag¬ 
netic field along the z-direction, which is perpendicular 
to the Zeeman field, with |a| being the strength of the 
spin-orbit coupled field. With this specific effective mag¬ 
netic field fl(k) by setting 5 = 0, it has been revealed 
that both the conventiona l 29 ' 30 ^ and anomalou s 31 ' 38 ' 39 
D’yakonov-Perel’ (DP)^ spin relaxations can be realized 
with (|fi z (k)|) > Cl and (|ST_ (k) |) <C Cl, respectively. For 
the conventional situation, in the strong (weak) scatter¬ 
ing limit when (|n(k)|)r^ <C 1 ((|fl(k)|)r^ > 1), the spin 
relaxation time (SRT) r s is inversely proportional (pro¬ 
portional) to the momentum scattering time r)t. (...) 
here denotes the ensemble average. For the anomalous 
situation ; 31 ' 38 ' 39 it has been found that by tuning the in¬ 
teratom interaction, the transverse spin relaxation can 
be divided into four regimes: the normal weak scatter¬ 
ing (r s oc t£), the anomalous DP-like (r s _1 oc r)*), the 
anomalous Elliott-Yafet (EY)-lik o 41 ' 42 (r s oc r^) and the 
normal strong scattering (t s -1 oc t^) regimes. Whereas 
the longitudinal spin relaxation can be divided into two: 





2 


i.e., the anomalous EY-like and the normal strong scat¬ 
tering regimes^ 

In contrast to the spin relaxation, the study for the 
spin diffusion in cold atoms with SOC has not yet been 
reported in the literature. However, the experimen¬ 
tal configuration realized by Brantut et al.,— for the 
“charge” diffusion of cold atoms can be applied to the 
spin diffusion straightforwardly. In their experiment, 
by adding the “barrier” laser in the middle of the cold 
atoms, the system is separated to the left and right parts. 
In the left part, the spin-polarized cold atoms can be 
prepared^ -28 ! 32 whereas in the right part, the system 
remains in the equilibrium state. Specifically, the atom 
densities in the left and right parts are prepared to be the 
same. With the SOC introduced to the right part ,— — & 
i.e., the spin diffusion region, by removing the “barrier” 
laser, this configuration can be used to study spin diffu¬ 
sion along one direction for the three dimensional (3D) 
Fermi gas with SOC. 

In the coordinate defined in Eq. ©. there are two 
specific configurations with the spin diffusions along the 
x- and y-directions, respectively. Accordingly, in the 
scattering-free situation, the k-dependent spin preces¬ 
sion frequencies in the spacial domain, i.e., the inhomoge¬ 
neous broadening in the steady-state spin diffusion 
along the x- and y-directions are determined by 

w x (k) = mCl(k)/k x = (riifl/k x ,0,ma S j , (2) 

u) y (k) = mtt(k)/k y = (^m£l/ky,0,mak x /ky S J , (3) 

respectively 52 i 53 Here, m is the atom mass. From Eq. © 
[Eq. ©], it can be seen that in contrast to the spin relax¬ 
ation in the time domain, in inhomogeneous broadening 
in the spacial domain, the original k-independent Zeeman 
field D becomes k-dependent, whereas the k-dependent 
spin-orbit coupled field mak x becomes k-independent 
(remains k-dependent). Hence, for the spin diffusion 
along the x-direction, the inhomogeneous broadening 
^(k) is similar to the one for spin relaxation in time 
domain fi(k) with one k-independent magnetic field per¬ 
pendicular to another k-dependent one. Accordingly, 
rich regimes may exit in the steady-state spin diffusion 
along the x-direction for both the spin polarizations per¬ 
pendicular and parallel to the Zeeman field. Whereas for 
the spin diffusion along the y-direction, from the differ¬ 
ent inhomogeneous broadenings in Eqs. © and m , its 
behavior should be very different from the one along the 
x-direction. 

In the present work, we investigate the steady-state 
spin diffusion for the 3D ultracold spin-orbit coupled 
40 K gas by the kinetic spin Bloch equation (KSBE) 
approach^ 4 both analytically and numerically. Four con¬ 
figurations, i.e., the spin diffusion along the x- and 
y-directions for the spin polarization P perpendicular 
(P || z, transverse configuration) and parallel (P || x, lon¬ 
gitudinal configuration) to the Zeeman field are studied. 
It is shown analytically that the behaviors of the steady- 
state spin diffusion for the four configurations are very 


different, which are determined by three characteristic 
lengths: the mean free path l T = kT^/m, the Zeeman os¬ 
cillation length Iq = k/(V3mCl) and the SOC oscillation 
length l a = l/(m|a|). The spin diffusion lengths for the 
spin diffusions in the four configurations are derived in 
the strong scattering regime, which are then extended to 
the weak scattering one. We further find that by dividing 
the system into different regimes, the complex analytical 
results can be reduced to extremely simple forms. It is 
revealed that by tuning the scattering, the system can 
be divided into five regimes: I, weak scattering regime 
(It > In,la)', II, Zeeman field-dominated moderate scat¬ 
tering regime (In <C l T la)', HI, SOC-dominated mod¬ 
erate scattering regime (l a -C l T "C In)', IV, relatively 
strong scattering regime (If -C l T "C ln,l a )', V, strong 
scattering regime (l T <C In, la, If)- Here, If represents the 
crossover length between the relatively strong and strong 
scattering regimes. In different regimes, the behaviors of 
the spacial evolution of the steady-state spin polarization 
are very rich, showing different dependencies on the scat¬ 
tering strength, Zeeman field and SOC strength. These 
dependencies are summarized in Tables [T] and eh for the 
spin diffusions along the x- and y-directions, respectively. 

The rich behaviors of the spin diffusions in differ¬ 
ent regimes are hard to be understood in the frame¬ 
work of the previous simple drift-diffusion model^^r— 
or the direct inhomogeneous broadening [Eqs. © and 
©] pictur o 44 : 51- — in the literature. In the simple drift- 
diffusion model, there are only two rather than five 
regimes: the strong scattering regime with l s oc l/(m|o:|) 
and the weak scattering regime with l s oc k^/rff/m. In 
the present work, it is found that the behaviors of the 
spin diffusions can be analyzed in the situation either 
with strong Zeeman and weak spin-orbit coupled fields 
(Regimes II and V) or weak Zeeman and strong spin- 
orbit coupled fields (Regimes III and IV). Accordingly, 
our previous inhomogeneous broadenings [Eqs. © and 
©] should be extended to the effective ones. It is found 
that when the spin polarization is parallel to the larger 
field between the Zeeman and spin-orbit coupled fields, 
the spin polarization cannot precess around the effective 
inhomogeneous broadening fields efficiently. In this situ¬ 
ation, the previous drift-diffusion model is applicable but 
r s (k) modified as follows ( modified drift-diffusion model). 
In the strong scattering regime, T s (k) remains the SRT in 
the conventional DP mechanism ; 29 : 30 ! 32 ! 44 whereas in the 
moderate scattering regime, T s (k) is replaced by the helix 
spin-flip rates determined in the helix space d 1 ' 39 When 
the spin polarization is perpendicular to the larger field 
between the Zeeman and spin-orbit coupled fields, the 
spin polarization can rotate around the effective inho¬ 
mogeneous broadening fields efficiently. Hence, the be¬ 
havior of the spin diffusion is determined by the effec¬ 
tive inhomogeneous broadenings together with the spin- 
conserving scattering (modified inhomogeneous broaden¬ 
ing picture). Based on the modified drift-diffusion model 
and modified inhomogeneous broadening picture, apart 
from Regime IV, all the features in different regimes can 
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be well obtained. 

Several anomalous features of the spin diffusion, which 
are in contrast to those obtained from both the sim¬ 
ple drift-diffusion model and the direct inhomogeneous 
broadening picture, are revealed. In the scattering 
strength dependence, it is found that when l a -C Iq, 
the longitudinal spin diffusion along the y-direction is 
robust against the scattering in a wide range including 
both the strong and weak scattering regimes. In the Zee- 
man field dependence, when the system is in Regime II, 
the longitudinal spin diffusion is enhanced by the Zeeman 
field. In the SOC strength dependence, we find that the 
spin diffusion length can be also enhanced by the SOC 
in Regime III. All these anomalous behaviors have been 
well understood from our modified drift-diffusion model 
and/or modified inhomogeneous broadening picture. 

This paper is organized as follows. In Sec. El we set 
up the model and KSBEs. In Sec. EH we show the an¬ 
alytical results for the transverse and longitudinal spin 
diffusion along the x- and y-directions. Different depen¬ 
dencies of the spin diffusion length on the mean free path, 
Zeeman oscillation length and SOC oscillation length are 
revealed. In Sec. hyi both the analytical and numer¬ 
ical calculations for the steady-state spin diffusion in 
3D isotropic speckle disorder are presented. Specifically, 
the disorder strength (Sec. II V Al) . Zeeman field strength 
(Sec. IIVBD and SOC strength (Sec. II V Cl) dependencies 
are discussed. We conclude and discuss in Sec. El 


II. MODEL AND KSBES 

With the 3D disordered speckle potential introduced 
to the spin diffusion region^ - — the Hamiltonian of 
the spin-orbit coupled ultracold atom, which is com¬ 
posed by the effective Hamiltonian Hq the disor¬ 
dered speckle potential U (r)^ ~ 37 i 43 and the interatom 
interaction H ln t, is written as 

H = H 0 + U(r)+H int . (4) 

The effective Hamiltonian consists of the kinetic energy 
of the atom and the SOC (h = 1), 

H 0 = k 2 / (2m) + «(k) • er/2, (5) 

with er being the vector composed of the Pauli matri¬ 
ces. The interaction Hamiltonian H ln t is approximated 
by the s-wave interatom scatteringI n our study, 
the scattering length is tuned to be zero by the Fesh- 
bach resonance/ ^ 37 ! 43 and hence Hint is absent in our 
following discussion. 

The KSBEs, derived via the nonequilibrium Green 
function method with the generalized Kadanoff-Baym 
AnsatZ ) 44 i 51 i 59 ~ — are utilized to study the spin diffusion 
in the ultracold Fermi gas: 

d t pk(r,t) = 3 t p k (r,t)|dif + d t pk(r,t)\ coh + <9 i( o k (r, t)| scat . 

( 6 ) 


In these equations, p k (r, t) represent the density matrices 
of atom with momentum k at position r = (x, y , z) and 
time t, in which the diagonal elements pk.aa describe the 
atom distribution functions and the off-diagonal elements 
Pk,a-<j represent the correlation between the spin-up and 
down states. 

For the quasi-one dimensional spin diffusion, the diffu¬ 
sion term is written as 

dtpk{r,i)\ A is = -(fcf/m)d c p k (r,t), (7) 

with ( = x or y for the spin diffusion along the x- or 
y-direction, respectively. The coherent term is given by 

^tPk(r,t)|coh = i[fl(k) • <r/2,p k (r,i)], (8) 

where [ , ] denotes the commutator. 

The scattering terms 9 t p k (r, i)| scat represent the atom- 
disorder scattering; 3 ' In our study, the effective Zeeman 
splitting energy and the SOC energy are much smaller 
than the Fermi energy^ Hence the atom-disorder scat¬ 
tering reads 

dtPkltL = 2tt^ |D k _ k / 2 <5(e k - e k ')(Pk' - Pk), (9) 

k' 

where 

|D q | 2 = J j dvdv'{[U(v)-Uo}[U{v')-Uo])e-^ r -^ 

= J y’drdr , C'(r-r , )e-' q (r - r ' ) = C q . (10) 

In Eq. (Hhd, Do is the average value of the disorder po¬ 
tential. For the 3D isotropic disordered speckle^ 

C q = 7T 3/2 V^cr|exp(-cr|( ? 2 /4), (11) 

where Vr is the potential amplitude and or denotes the 
radius of the auto-correlation function of the laser. 

In our model, with the same atom densities and hence 
the same chemical potentials in the left and right parts of 
the system, there is no “charge” diffusion in the system. 
The spin polarization at the boundary between the left 
and right parts of the system is approximately treated to 
be fixed. 


III. ANALYTICAL RESULTS 


In this section, we analytically study the steady-state 
spin diffusion along the x- and y-directions in cold atoms 
with the atom-disorder scattering [Eq. ©]. Both the sit¬ 
uations with spin polarization perpendicular (P||z) and 
parallel (P||x) to the Zeeman field are analyzed. 

In the steady state, the KSBEs are written as 


(k c /m)d c pk( r) + i Qa x / 2, p k (r) 
+ Wkk’ [pk(r) - Pk'( r )] = 0 , 


+ i 


ak x a z /2,pk(r) 

(12) 
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where Wkk' = 27rCk-k"5(£k — £k')- In the strong scat¬ 
tering regime with l T -C In, la [he., (|fi(k)|)Tfc -C 1], the 
steady-state diffusion lengths for different configurations 
are obtained and extended to the situation with moder¬ 
ate scattering strength (Zn -C l T <C l a or l a -C l T -C Zn)- 


The corresponding oscillation length 


In, 
la , 

\J2lrln/V 3 , 


when In <C l T <C l a 
when l a € I T € Zn 

when l T <C Z a &Z r <C Zn 


(17) 


A. Spin diffusion along the x-direction 


When the spin diffusion is along the x-direction, by 
taking the steady-state condition, the Legendre com¬ 
ponents of the azimuth-angle-averaged density matrix 
[Eq. (I A2[) ] with respect to the zenith angle 6 k are given 
by 


k d r lp‘ k 


i -i 


(l + 1 )p\ 


Z+l _ 


m dx V y/21 - 1 


+ i 


ak 

l~2 l 


i -1 
'k 


V2T+3 
, (i +1 )p! 


rf2 r- - 

*[-^0'*,PfcV2Z +1 


^+i n 
k 


VzT^l v/ 2T+3 j 


-^v / 2T+T = 0, 

Tfc f 


(13) 


in which the momentum relaxation time is denoted by 

C^cos 9) [l — Pi (cos 0)] sin Odd, (14) 

with Pi (cos 9) being the Legendre function. By keeping 
both the zeroth and first orders (Z = 0,1), the analyt¬ 
ical solution for the spin polarization is obtained from 
Eq. (fl3ll (refer to Annendix IA 11) . It is found that for both 
the spin polarization in the transverse (x-T) or longitu¬ 
dinal (x-L) configuration, the spin polarization is limited 
by one oscillation decay together with one single expo¬ 
nential decay, i.e., 

S| ~ A£ exp (—x/L x ) cos (x/l x ) + Be exp(— x/L x ), (15) 

with £ = t (transverse) or l (longitudinal). In Eq. (fT5l) . 
A ,£ and B^ are the amplitudes for the oscillation and 
single exponential decays, respectively, which are de¬ 
termined by the boundary condition; L x and L x are 
the decay lengths for the oscillation and single expo¬ 
nential decays, respectively; l x is the oscillation length 
for the oscillation decay. The integral forms for L x , L x 
and l x are complicated [Eqs. (1A9[) , (IA10I) and (I A11 ft in 
Appendix [All . However, when the system is in the 
strong (l T -C In, la) and moderate (Zn <C l T <C l a or 
l a <C l T <C In) scattering regimes, it is found that the 
analytical results can be reduced to simple forms, which 
can describe the behavior of the spin diffusion quite well 
(Sec. 13). 

Specifically, for the oscillation decay, 

{ 2Z T /v / 3, when In <C Z T <C l a 
V 2l T ln/(V3 l a ), when l a <C l T < In . (16) 

\j2l T ln/V3, when l T -C Z a &Z T < In 


mVk f 77 

Tk ’ 1 = ~ 2 ^~ J 0 


From Eq. m, one further notes that In and l a corre¬ 
spond to the spacial oscillation length due the Zeeman 
and spin-orbit coupled fields. Because of this, we refer 
to In and l a as Zeeman and SOC oscillation lengths, re¬ 
spectively. For the single exponential decay, the diffusion 
length reads 

! Zi-Za/ (V$ln), when In It la 

Irln/iVSla), when l a < l T < Zn ■ (18) 

l a , when l T <C Z a &Z r <C In 

From these simple dependencies of the spin diffusion 
length on the mean free path, the Zeeman oscillation 
length and the SOC oscillation length, the behavior of 
the steady-state spin polarization shown in Eq. m can 
be further simplified. It can be demonstrated that when 
L x « L x , Aj « B^. Specifically, only when L x ss L x , 
both the oscillation decay and single exponential decay 
are important, but with similar decay length. Other¬ 
wise, the spin polarization can be approximately reduced 
to one oscillation or single exponential decay, which de¬ 
pends on its amplitude in the spin polarization. In differ¬ 
ent regimes, the different behaviors for the steady-state 
spin polarization are analyzed as follows. In the Zeeman- 
field (SOC) dominated moderate scattering regime with 
Zn < It < l a (l a < It -C In), the condition L x ss L x s 
is never satisfied. Accordingly, when Zn -C Z r -C l a 
(l a <C l T ■C In), in the transverse/longitudinal situation, 
the steady-state spin polarization is approximately oscil¬ 
lation/single exponential (single exponential/oscillation) 
decay. In the strong scattering regime (l T <C l a ,ln), when 
L x ~ L x , one obtains 

Zr,* » Vsi 2 J(2l n ), (19) 

which is referred to as the crossover length between the 
relatively strong and strong scattering regimes. Accord¬ 
ingly, when l T » Z£ (Z T <C l c 7 : ), the steady-state 
spin polarization is approximated by a single exponen¬ 
tial/oscillation (oscillation/single exponential) decay in 
the transverse/longitudinal situation. 

Based on the above analysis, we summarize the be¬ 
haviors of the steady-state spin polarization and the spin 
diffusion lengths in Table Q] for the two specific config¬ 
urations x-T and x-L in different regimes. As shown 
in the table, the system is divided into five regimes: I, 
weak scattering regime (l T In, la)’, II, Zeeman field- 
dominated moderate scattering regime (In <C l T 
la)’, HI, SOC-dominated moderate scattering regime 
(l a C Z T C In); IV, relatively strong scattering regime 
(1% X <C Z T <C l a , In); V, strong scattering regime (l T -C 
la, In, It x)- I n different regimes, it can be seen that the 
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dependencies of the spin diffusion length on the scatter¬ 
ing strength, the Zeeman field and SOC strength are very 
rich. Specifically, in the Zeeman field- (SOC-) dominated 
moderate scattering regime with Iq -C l T -C l a (l a -C 
l T In), the longitudinal (longitudinal/transverse) spin 
diffusion is enhanced by the Zeeman field (SOC); whereas 
in the relatively strong (strong) scattering regime (l T <C 
la, In), the transverse (longitudinal) spin diffusion is de¬ 
termined by only the SOC oscillation length, but irrele¬ 
vant to the Zeeman field with l T 1° x (l T <C l° x )- The 
rich behaviors of the spin diffusions in different regimes 
are hard to be understood in the framework of the previ¬ 
ous simple drift-diffusion modeler— or the direct inho¬ 
mogeneous broadening [Eqs. (0 and 0] pictur ed 4 ’ 51- — 
in the literature. In the simple drift-diffusion model, 
there are only two rather than five regimes: the strong 
scattering regime with l s oc l/(ro|a|) and the weak scat¬ 


tering regime with l s oc ky/rk/m. In the direct inhomoge¬ 
neous broadening picture, from Eq. 0, the Zeeman field 
(SOC) can (cannot) provide the inhomogeneous broad¬ 
ening. Hence, it seems that the spin diffusion can be sup¬ 
pressed by the Zeeman field, but irrelevant to the SOC. 
Below, we extend our previous inhomogeneous broaden- 
ings [Eqs. 0 and 0] to the effective ones. We will show 
that based on the effective inhomogeneous broadening, 
apart from Regime IV, the above anomalous behaviors 
of the spin diffusion can be understood from the view 
point of the helix representation ; 31 ! 39 In the helix space, 
apart from the spin-conserving scattering [the fifth term 
in the left-hand side of Eq. (A6) in our previous worl—.], 
additional terms arise including the helix spin-flip scat¬ 
tering [the sixth term in the left-hand side of Eq. (A6) 
in Ref. [31] and helix coherence term [the last term in the 
left-hand side of Eq. (A6) in Ref. [3l| ^ 


TABLE I: Behaviors of the steady-state spin polarization in the spatial domain and corresponding spin diffusion lengths for 
configurations x-T and x-L in different regimes. 


Regime 

Condition 

Behavior and Lf. in x-T 

Behavior and Lf, in x-L 

I: weak scattering regime 

It Ifl i l a 

NA 

NA 

II: Zeeman field-dominated 

moderate scattering regime 

In < l T < l a 

oscillation decay; 

2 l T /V 3 

single exponential decay; 
lrla/(V 3 ln) 

III: SOC-dominated moderate 

scattering regime 

LCItC In 

single exponential decay; 

l T ln/(V 3 l a ) 

oscillation decay; 

V 2 l T ln / (V 3 l a ) 

IV: relatively strong scattering 

It « l a ,ln&dr > 1%, X 

single exponential decay; 

oscillation decay; 

regime 

(It < la < In) 

la. 

\J 2 l r ln/V 3 

V: strong scattering regime 

It < la,ln&d T < l C T, x 

oscillation decay; 

\J 2l T ln/V3 

single exponential decay; 

la 


It,x ~ V3ll/(2ln). 


When the system is in Regime II, the Zeeman field- 
dominated moderate scattering regime (In -C l T ■C la) or 
Regime V, the strong scattering regime (l T <C l a ,ln -C 
1% X ), the condition In -C l a is satisfied. Therefore, both 
the transverse and longitudinal spin diffusions can be 
understood in the limit with strong Zeeman and weak 
spin-orbit coupled fields. In this situation, the effective 
inhomogeneous broadening field is given by 

Wg ff (k) = (m/k x )\JVL 2 + a 2 k 2 x 

ss [mCl/k x + ma 2 k x / (2£l)] x', (20) 

with x' = --1-x + . akx / n -z being 

\/l + (ak x /ti ) 2 y/l + (ak x /Q ) 2 
nearly parallel to x. Under this effective field, the be¬ 
haviors of the spin precession in the spacial domain for 
the transverse and longitudinal spin diffusions can be ob¬ 
tained, as schematically shown in Fig. [1] In this figure, 


for the transverse spin diffusion, the spin polarization 
is perpendicular to the strong Zeeman field and hence 
Wg ff (k) approximately. During the scattering, the spin 
vectors rotate around the effective inhomogeneous broad¬ 
ening field CLV ff (k) fast. Moreover, the scattering can 
also influence the spin diffusion. In the helix space, as 
mentioned above, apart from the original spin-conserving 
scattering, there exist the helix spin-flip scattering and 
helix coherence processes. However, with the helix spin- 
flip rate a 2 k 2 /(fl 2 Tk) -C 1 /rfc and helix coherence rate 
ak/fflrff) <C 1/Vfc when H ak , both the helix spin- 

flip scattering and helix coherence can be neglected. In 
this situation, the effective inhomogeneous broadening 
together with the spin-conserving scattering determines 
the behavior of the spin diffusion. We refer to this pic¬ 
ture as the modified inhomogenous broadening picture. 
For the longitudinal situation, the spin polarization is 
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FIG. 1: (Color online) Schematic for the spin precession 
around the effective inhomogeneous broadening field cuj ff (k) 
[Eq. (1201) 1 in the transverse (Sj) and longitudinal (S£) spin 
diffusions. With the strong Zeeman and weak spin-orbit cou¬ 
pled fields, for the transverse (longitudinal) situation, the spin 
vector Sj (S£) is perpendicular (parallel) to tUg ff (k) approx¬ 
imately. Therefore, the inhomogeneous broadening held can 
(cannot) cause efficient spin precession in the transverse (lon¬ 
gitudinal) spin diffusion. 


nearly parallel to u>g ff (k), and hence the effective inho¬ 
mogeneous broadening cannot cause the spin precession 
effectively. In this situation, the spin diffusion can be 
understood from the drift-diffusion mo del^ - — modified 
as follows. The diffusion length l s = yjD t s (k) in which 
D = VpTk/3 is the diffusion coefficient with vf being 
the Fermi velocity. The SRT is analyzed in the helix 
space. First of all, for the longitudinal situation here, 
the helix coherence term has no contribution to the spin 
relaxation. In this situation, there exit two channels 
influencing the spin relaxation: (i), the effective inho- 
mogenous broadening together with the spin-conserving 
scattering; (ii), the helix spin-flip scattering ] 31 ' 39 In the 
strong scattering regime, both channels (i) and (ii) are 
important for the spin relaxation, in which t s (k) remains 
the SRT in the conventional DP mechanism j 29 i 30 i 32 i 44 In 
the moderate scattering regime, channel (ii) is dominant 
for the spin relaxation, and hence t s (k) is replaced by the 
helix spin-flip rates determined in the helix space . 31 ' 39 
We refer to the above pictures as modified drift-diffusion 
model. Accordingly, in the moderate scattering regime, 
the dependence of the helix spin-flip rate on the scat¬ 
tering strength, the Zeeman field and SOC strength dis¬ 
closed in our previous work s 31 ' 38 ' 39 can also influence the 
spin diffusion. 

Specifically, in the Zeeman field-dominated moderate 
scattering regime (In < 1 T C l a ), i.e., Regime II, the 
Zeeman oscillation length is the shortest length scale 
in the spin diffusion, which determines the behavior of 
spin precession in the spacial domain during the scat¬ 
tering. When the spin polarization is perpendicular to 
the Zeeman field (transverse configuration), the spin vec¬ 
tor approximately precesses around mfl/Arx', leading to 


the spacial oscillations with the period proportional to 
(\k x \ )/(mSl) K ~ in Eq. iHTjl ]. Moreover, due to the 
fast spacial oscillations with the strong Zeeman field, the 
spin memory is lost during one spin-conserving scatter¬ 
ing, with the diffusion length being approximately the 
mean free path (Uf « 2l T /y/3 in Table Q]). When the spin 
polarization is parallel to the Zeeman field (longitudinal 
configuration), the effective inhomogeneous broadening 
Wp ff (k) cannot cause spin precession efficiently and the 
steady-state spin polarization decays without any oscil¬ 
lation. The spin diffusion can be understood from the 
modified drift-diffusion modeller— Specifically, the he¬ 
lix spin-flip rate is calculated to be a 2 k 2 /(2£l 2 t^) in the 
moderate scattering situation . 31 " 38 " 39 Accordingly, the 
spin diffusion length in the modified drift-diffusion model 
is given by l s ~ \f2l T l a /(V^>ln)i which is consistent with 
our model shown as L x l « l T la/(V3In) hi TableQ] Specif¬ 
ically, one notes that due to the suppression of the spin 
relaxation by the Zeeman field, the longitudinal spin dif¬ 
fusion length is enhanced by the Zeeman field. 

In Regime V, the strong scattering regime ( l T -C 
Io.iIq.iIt x)i we consider a limit situation with the Zee- 
man field much stronger than the spin-orbit coupled one 
(l T l a , <C l c TX ). The effective inhomogeneous broad¬ 
ening is given by Eq. For the transverse spin dif¬ 

fusion, because the inhomogeneous broadening given by 
the Zeeman field is dominant, the spin diffusion length is 
suppressed by the Zeeman field, but less influenced by the 
SOC. Moreover, during the diffusion, the atoms experi¬ 
ence several spin-conserving scatterings, which suppress 
the spin diffusion. Accordingly, we obtain a reasonable 

picture to understand ^ J 2Mfi/v / 3 in Table in For 

the longitudinal spin diffusion, the inhomogeneous broad¬ 
ening cannot cause the spin precession efficiently, with 
the steady-state spin polarization showing single expo¬ 
nential decay. Hence, the modified drift-diffusion model 
can be used— - — With the SRT in the strong scattering 
regime r s (k) a: 2/(a 2 fc 2 7fc) ) 31 i 38 i 39 one obtains the spin 
diffusion length l s ss \/2l a /V3 (if. ~ l a in Table Q]). 
Therefore, the longitudinal spin diffusion length depends 
only on the SOC oscillation length] 31 ' 45- — 

When the system lies in Regime III, the SOC- 
dominated moderate scattering regime (l a <C l T In) 
and Regime IV, the relatively strong scattering regime 
(1% x -C l T lajn)i the spin-orbit coupled field is much 
stronger than the Zeeman one. In this situation, the ef¬ 
fective inhomogeneous broadening field reads 


w'e ff (k) = (m/k x )^a 2 kl + D 2 z' 

~ \ma + mQ 2 /(2ak 2 )\ z', (21) 


where z ' 


] , n/(akx) is 

■\/l + [0/(a/c x )] 2 \/l V [0/ (ak x )} 2 


parallel to z approximately. The analysis is similar to 
the situation with strong Zeeman and weak spin-orbit 
coupled fields. One obtains that for the transverse (lon¬ 
gitudinal) spin diffusion, the spin polarization is nearly 
parallel (perpendicular) to u/g ff (k), which cannot (can) 
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cause efficient spin precession. These pictures are sum¬ 
marized in Fig. [2] Therefore, the behavior of the trans- 



FIG. 2: (Color online) Schematic for the spin precession 
around the effective inhomogeneous broadening field ui 'Sr(k) 
[Eq. (HUll l in the transverse (Sj) and longitudinal (S£) spin 
diffusions. With the weak Zeeman and strong spin-orbit cou¬ 
pled fields, for the transverse (longitudinal) situation, the spin 
vector Sj (S£) is parallel (perpendicular) to u/g ff (k) approx¬ 
imately. Therefore, the inhomogeneous broadening field can¬ 
not (can) cause efficient spin precession in the transverse (lon¬ 
gitudinal) spin diffusion. 

verse spin diffusion can be understood from the modified 
drift-diffusion modeljs* — whereas the longitudinal spin 
diffusion can be analyzed from the features of the preces¬ 
sion of the spin vectors around u/g ff (k). 

Specifically, in Regime III, i.e., the SOC-dominated 
moderate scattering regime (l a <C l T ^ In), in the trans¬ 
verse configuration, the modified drift-diffusion model is 
used to understand the spin diffusion^ - — In the he¬ 
lix representation, the helix spin-flip rate is proportional 
to ft, 2 /(a 2 k 2 Tk) approximately 31 ' 38 ! 39 The corresponding 
spin diffusion length is proportional to l T ln/la, which 
is consistent with Lf, ~ l T ln/(V^l a ) in Table Q] Con¬ 
sequently, one observes that the spin relaxation is sup¬ 
pressed by the spin-orbit coupled field, and the transverse 
spin diffusion length is enhanced by the SOC. In the lon¬ 
gitudinal configuration, with m\a\ mti 2 / {2\a\k 2 ) , the 
spin polarization evolves with oscillations in the spacial 
domain, whose oscillation length is l x ~ l a [Eq. (fTTl) ]. 
Furthermore, only mU 2 / (2ak 2 )z’ in u/g ff (k) can cause 
the inhomogeneous broadening. Therefore, it can be ex¬ 
pected that the longitudinal spin diffusion length is in¬ 
versely proportional (proportional) to the Zeeman field 
(SOC). Moreover, the spin-conserving scattering can sup¬ 
press the spin diffusion in the modified inhomogeneous 
broadening picture. These analysis are consistent with 
14 « y/2l r l n /(V3l a ) in Table U 

In Regime IV, the relatively strong scattering regime 
(1% x <C l T -C l a , In), for the transverse spin diffu¬ 
sion, from the modified drift-diffusion model^~— l s ~ 
V2L/V3 (L x t « l a in Table |U) . For the longitudinal spin 
diffusion, one expects that the modified inhomogeneous 


broadening picture can be applied. However, this pic¬ 
ture fails to explain the behavior of the longitudinal spin 
diffusion. It can be seen from Eq. (ED that both the 
spin-orbit coupled and Zeeman fields provide the inho¬ 
mogeneous broadening [mCl 2 / (2ak 2 )z' in Eq. ED]- By 
further considering that the spin diffusion is suppressed 
by the spin-conserving scattering, it is obtained that the 
spin diffusion length is proportional to the SOC strength, 
but inversely proportional to the Zeeman field strength 
and scattering. However, in Table Q] the longitudinal 

spin diffusion length L x l ss \J 21 t Iq, / y/H is irrelevant to 
the SOC strength. One further notices that Regime IV 
lies in the crossover region between the moderate and 
strong scattering regimes. When the scattering is rela¬ 
tively strong, the shortest length scale in the spin diffu¬ 
sion is the mean free path. However, there still exists 
strong competition between the effective inhomogeneous 
broadening and scattering, which makes the behavior of 
the spin diffusion complicated ! 50 ! 53 ! 54 

Finally, we address the behavior of the spin diffusion 
along the x-direction in the limit situation where the Zee- 
man field is zero. When = 0, In is infinite. The corre¬ 
sponding steady-state spin polarization shows very differ¬ 
ent behaviors compared to the situation with finite Zee- 
man field. In this situation, the transverse (longitudinal) 
spin diffusion length is infinite because (k) cannot 
cause inhomogeneous broadening^ Specifically, for the 
longitudinal situation, the spin polarization is perpendic¬ 
ular to the spin-orbit coupled field, spin helix establishes 
with the oscillation length being l a -—^ 


B. Spin diffusion along the y-direction 


When the spin diffusion is along the y-direction, 
by taking the steady-state condition, the Fourier com¬ 
ponents of the azimuth-angle-averaged density matrix 
[Eq. (I A2[) ] with respect to the zenith angle 6 k are given 

by 


hk d 
2im dy 

r n 




ak _ ;+1 
l~ 4 -Vz’Pk 


~ i -1 

Pk 


+ i 


L 2 


Pk 


+ ^ = 0 , 
T~k,l 


( 22 ) 


in which the momentum relaxation time is denoted by 


C(cos0)[l — cos(Z0)] sin Odd. (23) 

One notes that by choosing the Legendre and Fourier 
expansions of the density matrix for the spin diffusions 
along the x- and y-directions, the definitions for the mo¬ 
mentum scattering time Tk,i [Eq. (fldl) ] and f k,i [Eq. ED ] 
are different. However, when 1 = 1, the two definitions 
are the same. By keeping both the zeroth and first orders 
(l = 0,1), the analytical solutions for both the spin po¬ 
larization perpendicular (y-T) and parallel (y-L) to the 
Zeeman field are obtained (refer to Appendix IA 21) . 


--1 


iy/k 


2n 






















When the spin polarization is perpendicular to the Zee- 
man field, in the moderate and strong scattering regimes, 
similar to the spin diffusion along the x-direction (Ta¬ 
ble |U), the behaviors of the steady-state spin polarization 
are different with different scattering strengths. When 
the scattering is relatively weak, which satisfies l T > 1° 

with , the steady-state spin 

v3fo 3 /q/ 

polarization for is limited by the bi-exponential decay, 
i.e., 


= P + exp(-y/L% + ) + P exp (~y/L% ), (24) 

with Ly ± being the diffusion length. It is further demon¬ 
strated that when Ly + <C Lji ~, P + <C P~ and hence the 
spin polarization [Eq. EH) ] reduces to a single exponen¬ 
tial decay with the decay length being L^~. Specifically, 
in the moderate scattering regime with l a <C l T •C In 
(In < l T <C l a ), the condition l T > is natu¬ 
rally (never) satisfied; in the strong scattering regime 
(l^ y < l T <C l a , In), it can be obtained that l a -C In and 
hence l% y « 4Zj(;/(v / 3Zn). Accordingly, when l T l^ y , 
the diffusion length for S V T is written as 

{ NA, when In -C l T <C l a 

V3l T ln/(2l a ), when l a <C l T < In 
V3l T ln/(2l a ), when l c T y < l T < l a , l n 

, (25) 

When the scattering is relatively strong, which satisfies 
l T < l%, y , the transverse spin polarization in the steady 
state is determined by the oscillation decay, i.e., 

Sy = P 0 exp (~y/L v T ) cos (y/l v T ). (26) 


Here, the decay length Lj. and oscillation length Zj. can 
be written as 


L v t 


\ V2l T , 

1 NA ’ 
[\/ VSItIq, 


when In <C l T l a 

when l a <C l T In 

when l T < l a , In, It, y 


(27) 


and 


l y T 


VV2ln, 

NA, 


when In <C l T <C l a 
when l a -C l T -C In , 
when l T < l a , In, It, y 


(28) 


respectively. 

When the spin polarization is parallel to the Zeeman 
field, it is found that the steady-state spin polarization 
is limited by the single exponential decay, i.e., 


S y L = P 0 exp(-y/L v L ). (29) 


Here, L v l = las/l^K 3Z^) + 1 is the decay length for S y L . 
In the moderate and strong scattering regimes, 


I Irla/ ()? when Iq It la 
L v l « < l a , when l a <C l T <C In ■ (30) 

[ l a , when l T <C Zn&Z r <C l a 

Based on the above results, in different regimes, the 
behaviors of the steady-state spin polarization and dif¬ 
fusion lengths are summarized in Table |TT1 for the two 
specific configurations y-T and y-L. 


TABLE II: Behaviors of the steady-state spin polarization in the spacial domain and corresponding spin diffusion lengths for 
configurations y-T and y-L in different regimes. 


Regime 

Condition 

Behavior and L ^ in y-T 

Behavior and L v l in y-L 

I: weak scattering regime 

It Id) la. 

NA 

NA 

II: Zeeman field-dominated 

moderate scattering regime 

IQ. <C It < la 

oscillation decay; 
y/2 l T 

single exponential decay; 
l T la/(VSln) 

III: SOC-dominated moderate 

scattering regime 

la < It < In 

single exponential decay; 

V3l T ln/(2l a ) 

single exponential decay; 

la 

IV: relatively strong scattering 
regime 

It < la,ln&d T > lr,y 

(It < la < In) 

single exponential decay; 

V3l T ln/(2l a ) 

single exponential decay; 

la 

V: strong scattering regime 

It < la,ln&d T < l C T,y 

oscillation decay; 

\J \[Zl T ln 

single exponential decay; 

la 


lr,y ~ 41„/(v / 31n). 


From Tabic ITI1 it can be seen that the spin diffusion along the y-direction is divided into similar five regimes as the 
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spin diffusion along the x-direction. The anomalous be¬ 
haviors for the spin diffusion along the x-direction also 
exist here. Nevertheless, new features arise in the spin 
diffusion along the y-direction. It is shown in Table [TT] 
that in Regimes III, IV and V, the longitudinal spin dif¬ 
fusion length is only determined by the SOC oscillation 
length, but irrelevant to the scattering. This robustness 
to the scattering for the spin diffusion in a wide range is 
further revealed in the scattering dependence of the spin 
diffusion (Sec. II V A[) . Below, it is found that based on 
the modified drift-diffusion model and modified inhomo¬ 
geneous broadening picture, apart from Regime IV, all 
the features in different regimes can be obtained. 

We first analyze Regime II, the Zeeman field- 
dominated moderate scattering regime (In <C l T "C la) 
and Regime V, the strong scattering regime (l T <C 
l a ,ln -C lr, y )- In these two regimes, with strong Zeeman 
and weak spin-orbit coupled fields, the effective inhomo¬ 
geneous broadening field is written as 

w eff( k ) = (m/k v )y/n 2 + a 2 k 2 x' 

ft [mfl/ky + ma 2 kl/(2flk y )]x', (31) 

with x' nearly parallel to x. In Regime II, i.e., the Zee- 
man field-dominated moderate scattering regime (In <C 
It "C la), when the spin polarization is perpendicular to 
the Zeeman held (transverse configuration), in the spacial 
domain, the spin vectors precess approximately around 
( mCl/ky )x', which causes spacial oscillations whose pe¬ 
riod is proportional to (\k y \)/(mCl) [Ij, ft \J3/2ln in 
Eq. ||2HD]. Moreover, due to the fast spacial oscillations, 
the spin memory is lost during one scattering, and hence 
the spin diffusion length is the mean free path approx¬ 
imately (Lj, ft a/2 l T in Table [Til) - When the spin po¬ 
larization is parallel to the Zeeman held (longitudinal 
configuration), the effective inhomogeneous broadening 
held a>g ff (k) cannot cause spin precession efficiently and 
the steady-state spin polarization decays without any os¬ 
cillation. From the modified drift-diffusion modeler— 
as calculated in Sec. MM Is » V2l T l a /(V3ln) [L v l ft 
l r la/(V3ln) in Table ITT). 

In Regime V, i.e., the strong scattering regime (l T -C 
la, In, It, y ), we analyze the limit situation with strong 
Zeeman and weak spin-orbit coupled helds (l T -C l a , In •C 
lr, y )- For the transverse spin diffusion, the inhomoge¬ 
neous broadening is dominantly determined by the Zee- 
man held [Eq. (13111 ]. Hence, the transverse spin diffu¬ 
sion length is suppressed by the Zeeman held, but less 
influenced by the SOC. Furthermore, during the diffu¬ 
sion, the spin-conserving scattering suppresses the spin 
diffusion (L v l ft \J\Z3l T ln in Table [TT1). For the longitudi¬ 
nal spin diffusion, the inhomogeneous broadening cannot 
cause the spin precession efficiently (single exponential 
decay). According to the modihed drift-diffusion model, 
the spin diffusion length depends only on the SOC oscil¬ 
lation length (Lj, ft l a in Table Q]). 

We then analyze Regime III, SOC-dominated moder¬ 
ate scattering regime (l a <C l T <C In), and Regime IV, 


relatively strong scattering regime (l( : < 1 T C l a ,ln)- 

With weak Zeeman and strong spin-orbit coupled helds, 
the effective inhomogeneous broadening held reads 

^'eff( k ) = (m/ky)^a 2 k%+n 2 z.' 

ft [mak x /ky + mfl 2 /(2ak x k y )\z r , (32) 

where z ' is parallel to z approximately. In the SOC- 
dominated moderate scattering regime (Regime III with 
la "C It <C In), in the transverse configuration, the spin 
polarization is nearly parallel to u/g ff (k). Therefore, 
the effective inhomogeneous broadening cannot cause 
the spin precession effectively (single exponential decay). 
From the modihed drift-diffusion model, the diffusion 
length is proportional to InW/la, which is consistent with 
Ly ft V3lnl T /(2l a ) in Table El In the longitudinal con¬ 
figuration, the steady-state spin polarization is perpen¬ 
dicular to c>/g ff (k) approximately. One notes that in 
uj'I s ( k), the spin-orbit coupled held (mak x /k y )z' pro¬ 
vides the dominant inhomogeneous broadening. More¬ 
over, this inhomogeneous broadening not only depends 
on k y but also k x , which can be more efficient than the 
one depending only on k x or k y . Due to this efficient 
inhomogeneous broadening, the steady-state spin polar¬ 
ization decays due to the interference without oscillation. 
Moreover, the spin memory can be lost in the scale of 
SOC oscillation length, which cannot persist in the mean 
free path (L v l ft l a in Table |TT|- This is different from the 
transverse spin diffusion in the Zeeman hled-dominated 
moderate scattering regime (In <C l T l a ). 

In the relatively strong scattering regime (Regime IV 
with l^ y < ! r < l a , In)- From Table El one observes 
that the behaviors for both the transverse and longitu¬ 
dinal spin diffusions are similar to the ones in Regime 
III, i.e., the SOC-dominated moderate scattering regime 
(l a <C l T <C In)- We address that this behavior is hard 
to be understood from both the modified drift-diffusion 
model and modified inhomogeneous broadening picture. 
For the transverse spin diffusion, the spin polarization is 
nearly parallel to (mak x /k y )i' , and hence the spin diffu¬ 
sion length is \f2l a /V?> from the modihed drift-diffusion 
model. For the longitudinal spin diffusion, the spin po¬ 
larization is perpendicular to the inhomogeneous broad¬ 
ening Held (mak x /ky)i' approximately, and hence the 
SOC can suppress the spin diffusion. Moreover, the spin- 
conserving scattering can suppress the longitudinal spin 
diffusion. From this analysis, the longitudinal spin dif¬ 
fusion length is suppressed by the SOC strength and 
scattering, but irrelevant to the Zeeman Held (modified 
inhomogeneous broadening picture). However, the pic¬ 
tures above fail to explain the transverse and longitu¬ 
dinal spin diffusion lengths for Regime IV in Table El 
with ft \Z3l T ln/(2la) and L v l ft l a - This is because 
for the transverse situation, although the modified drift 
diffusion model can explain the spin diffusion along the 
x-direction, it is too rough to consider the anisotropy be¬ 
tween the diffusions along the x- and y-directions in the 
relatively strong scattering regime. This was first pointed 
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out by Zhang and Wu in the study of the spin diffusion in 
graphene. 50 - For the longitudinal spin diffusion, when the 
scattering is strong, the shortest length scale in the spin 
diffusion is the mean free path. In this situation, there 
exists strong competition between the effective inhomo¬ 
geneous broadening and scattering, which makes the be¬ 
havior of the spin diffusion complicated j 50 i 53 i 54 

Finally, we emphasize that from TableCIl for the trans¬ 
verse spin diffusion along the y-direction, in Regimes 
III, IV and V, the spin diffusion lengths are irrelevant 
to the scattering. Therefore, with weak Zeeman and 
strong spin-orbit coupled fields (l a -C In), a specific sit¬ 
uation can be realized where the spin diffusion length is 
robust against the scattering except with the extremely 
weak scattering. It is noted that in the strong scattering 
regime, this feature was predicted in the simple drift- 
diffusion model^ - — and was also revealed in graphene 
by the KSBE approach^ In this work, we have further 
extended it into the weak scattering regime. 

IV. NUMERICAL RESULTS 

In the numerical calculation, the KSBEs are solved by 
employing the double-side boundary conditions, 53 

p k (C = 0, t) = '^ kt y f ki + f kf 2 cr • n, k c >0 
Pk(C = L,t) = /£, k ( < 0 

(33) 

where / ktT = {exp[(£ k — Ha) / (ksT)] + 1} 1 is the Fermi 
distribution function at temperature T, with fi -^ stand¬ 
ing for the chemical potentials determined by the atom 
density n Q =^ k Tr[p k ] and the spin polarization P( 0) in 
the left part of the system; n denotes the spin polariza¬ 
tion direction; / k is the Fermi distribution at equilibrium. 
When the system evolves to the steady state, one obtains 
the diffusion length from the spatial evolution of the spin 
polarization P(£) = )T k Tr[p k (£)<r • n ]/n a . 

Within the experimental feasibility by referring to the 
experiment by Wang et al the parameters are chosen 
as follows. The lowest two magnetic sublevels 19/2, 9/2) 
and |9/2, 7/2) are coupled by a pair of Raman beams 
with wavelength A = 773 nm and the frequency dif¬ 
ference uj/(2tt) = 10.27 MHz. The Raman detuning 
5 = oj z — ui is set to be zero by choosing the Zeeman shift 
u] z /(2n) = 10.27 MHz. The recoil momentum and en¬ 
ergy are set to be k r = fco/10 with fc 0 = 27r/A, and hence 
E r = k // (2m) = 27r x 83.4 Hz. In our study, the SOC 
strength varies from 0.5ao to 6ao with «o = —2 k r /m. 
The strengths of the Zeeman field H vary from !0E r to 
450.EV.. Furthermore, the Fermi momentum is set to be 
&f = 30 k r 2& It is noted that with these parameters, the 
condition that the Zeeman and SOC energies are much 
smaller than the Fermi energy is satisfied. 

Moreover, the temperature is set to T = 0.3Tp 
with Tp being the corresponding Fermi temperature^ 
With these parameters, the thermal deBroglie wave¬ 
length AdB = h/\J2-KmkB r P ~ 0.26 /mi. For the 3D 


isotropic speckle disorder, Vn/ks = 1250 nK and ctr = 
0.27 /mn i 34 i 35 With these disorder parameters, the mean 
free path l T ss 5 /mi. In our study, the strength of the dis¬ 
order strength V is tuned by the laser ,M ~ 37 i 43 One notes 
that when ( V/Vr ) 2 > 20, l T < A^b- According to the 
Ioffe-Regel condition for the Anderson localization^ the 
Anderson localization may become relevant. Neverthe¬ 
less, in our study, to compare with the analytical results 
in different regimes revealed in Sec. cm the numerical 
calculations are extended to ( V/Vr ) 2 20. 


A. Scattering strength dependence 

In this part, we study the scattering strength depen¬ 
dence of the steady-state spin diffusion of spin-orbit cou¬ 
pled 40 K gas in the 3D isotropic speckle disorder. The 
SOC strength is set to be ao and the spin polarization is 
chosen to be P = 20%. For the spin diffusion along the 
x-direction (y-direction), both the transverse and longi¬ 
tudinal spin diffusion lengths Lf. and (Lj, and L v l ) 
are shown in Figs. [3]) a) and (b) [Figs. (3])c) and (d)]. 

We first analyze the spin diffusion along the x-direction 
[Figs. [3D a) and (b)[. For the transverse spin diffusion, it 
can be seen from Fig. [3][a) that no matter the Zeeman 
field is strong with f l = 450.E r (the blue dashed curve 
with squares) or weak with fl = 10E r (the red solid 
curve with circles), the transverse spin diffusion length 
Lf. decreases with the increase of the disorder strength. 
The underlying physics can be understood as follows. We 
first calculate the characteristic lengths for the system de¬ 
fined in Sec. cm the SOC oscillation length l a ss 0.6 /mi; 
the Zeeman oscillation length In ~ 0.1 /im (4 /mi) 
for O = 450 E r (10 E r )\ when 0 = 450FA (10 E r ), the 
crossover length between the relatively strong and strong 
scattering regimes 1% X ~ \/3Z^/(2Zn) ss 3 /mi (0.08 /mi), 
i.e., (V/Vr)c ~ 2 [(V/Vr)c ~ 65]. Furthermore, it is cal¬ 
culated that when l T ss l a , (V/Vr) 2 ~ 8; when l T « In, 
(V/V R ) 2 « 50 [(V/Vr) 2 w 1] for fi = 450£; r (10 E r ). Ac¬ 
cordingly, with the increase of the disorder strength, the 
system experiences several regimes. For Q = 450£V., the 
regimes are approximately divided into 

fl: l T > In, L, when (V/V R ) 2 < 8 

\ II: In < It < len when 8 < (V/V R ) 2 < 50 , (34) 

[v : l T < ln,loc,lr,x when (^/Vr) 2 > 50 

which are labelled by the blue Roman numbers with the 
boundaries indicated by the blue crosses at the lower 
frame of Fig. [3][b). Therefore, when (V/Vr) 2 < 8, the 
system is in Regime I, and our calculation shows that 
the transverse diffusion length decreases with the in¬ 
crease of the disorder strength. When 8 < (V/Vr) 2 < 
50 [(V/Vr) 2 50], the system lies in Regime II 

(V), and from Table Cl one comes to Lf. ~ 2l T /y/3 

(L?/ ss \j2l T ln/V$) with the steady-state spin polariza¬ 
tion showing oscillation decay. Hence, the transverse spin 
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diffusion length decreases with the increase of the disor¬ 
der strength. One notes that the corresponding results 
calculated from the analytical formula Eq. (lAlOl) (the 
green dot-dashed curve) agree with the numerical ones 
in Regimes II and V in Fig. [3ta). For 17 = 10 E r , the 
regimes are approximately divided into 

fl: h Z In, l a , when {V/V R ) 2 < 1 

I HI :l a <lr< In, When 1 < ( V/V R ) 2 < 8 
| IV : 1% <l T < l n , L, when 8 < ( V/V R ) 2 < 65 ’ 
[V: l T < In, l a , 1% x , when {V/Vr) 2 > 65 

(35) 

which are shown by the red Roman numbers with the 


boundaries indicated by the red crosses at the upper 
frame of Fig. [31(a). Specifically, in Regime I when 
( V/Vr ) 2 < 1, it is shown that the transverse spin dif¬ 
fusion length decreases with the increase of the scatter¬ 
ing strength. In Regime III (V) with 1 < ( V/Vr ) 2 < 
8 [(V/Vr) 2 65], from Table [H it is obtained that 

Lf, ft l T ln/(V3la) {Lj, ft \J 2 l r ln/V3) with the steady- 
state spin polarization being single exponential (oscilla¬ 
tion) decay. Hence, the transverse spin diffusion length 
decreases with the increase of disorder strength. One 
further notices that in Fig. [3](a), the result calculated 
from the analytical formula Eq. (IA10I) (the orange dashed 
curve) agrees with the numerical one in Regime V. 




(V/V R ) 2 I V/Vr) 2 

FIG. 3: (Color online) Scattering strength dependence of the steady-state spin diffusion of spin-orbit coupled 40 K gas with the 
3D isotropic speckle disorder. The SOC strength a = «o and the spin polarization P = 20%. For the spin diffusion along 
the x-direction (y-direction), the transverse and longitudinal spin diffusion lengths L f* and L £ ( L ^ and L^) are shown in 
Fig s - EH a ) and (b) [Figs.EUc) and (d)], respectively. Situations with strong (17 = 450i? r ) and weak (17 = 10 E r ) Zeeman fields are 
calculated both analytically and numerically. The blue (red) crosses on the frames indicate the boundaries between different 
regimes represented by the blue (red) Roman numbers at the lower (upper) frame when 17 = 450i? r (107%). It is shown that 
the analytical calculations agree with the numerical ones in the relatively strong/strong scattering regime and crossover region 
between the relatively strong and moderate scattering regimes. 


For the longitudinal spin diffusion along the x- 


direction, it can be seen from Fig. [31(b) that no mat- 
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ter the Zeeman field is strong (450 E r , blue dashed curve 
with squares) or weak (10E)., red solid curve with cir¬ 
cles), with the increase of the disorder strength, the spin 
diffusion length decreases first and then becomes insen¬ 
sitive to the scattering. This can be understood as fol¬ 
lows. With the increase of the disorder strength, the 
corresponding regimes can also be divided according to 
Eq. (1M1) [Eq. (1351) ] when Q = 450 E r (10 E r ). Specifi¬ 
cally, when ( (V/Vr ) 2 < 8 [( V/Vr) 2 < 1] for = 450 E r 
(10 E r ), the system is in Regime I, with the longitudi¬ 
nal spin diffusion suppressed by the scattering. When 
8 < {V/V R Y < 50 [1 < (V/V R ) 2 < 8] for fi = 450 E r 
(10E r ), the system lies in Regime II (III), and one comes 
to L x l ft l T l a /(V3ln) [L x l ~ y/2l T ln/(y/3l a )] with the 
steady-state spin polarization showing single exponential 
(oscillation) decay. When (V/Vr) 2 50 [(V/Vr) 2 65] 

for H = 450 E r (10If r ), the system lies in Regime V, it is 
obtained that ft l a (single exponential decay). There¬ 
fore, the longitudinal spin diffusion is suppressed first 
and then becomes insensitive to the scattering with the 
increase of the disorder strength. Also in Fig. [3])b) , for 
H = 450 E r , it is shown that the results calculated from 
the analytical formula Eq. (IA9I) (the green dot-dashed 
curve) agrees with the numerical one in both Regimes 
II and V; for H = 10 E r , the analytical results (the blue 
dashed curve) agree with the numerical ones in Regime 
V. 

We then turn to the spin diffusion along the y-direction 
[Figs- EKc) and (d)]. For the transverse spin diffusion, it 
is shown in Fig. H c) that for both the strong (450 E r , 
blue dashed curve with squares) and weak (10 E r , red 
solid curve with circles) Zeeman fields, the spin diffu¬ 
sion length decreases with the increase of the disorder 
strength. It is calculated that for f 1 = 450 E r (10 E r ), 
« 8.3 /im (0.2 /im), i.e., (V/Vr) 2 c ft 0.6 \(V/Vr) 2 c ft 
Accordingly, when the Zeeman field is strong (fl = 
450 E r ), one can divide the regimes according to Eq. m- 
Specifically, when (V/Vr) 2 < 8, the transverse spin dif¬ 
fusion is calculated to be suppressed by the scattering. 
When 8 < (V/Vr) 2 < 50 [(V/Vr) 2 50], one obtains 

from Table HI! that L j, ft y/2l T (L j, ft \J\Z31qIt) with the 
steady-state spin polarization showing oscillation decay. 
Therefore, with the increase of the scattering strength, 
the transverse spin diffusion is suppressed. When the 
Zeeman field is weak (10-EV), the regimes are approxi¬ 
mately divided into 



f I: lr>lu,L, when (V/Vr ) 2 < 1 

I HI :l a <l T <ln, when 1 < (V/V R ) 2 < 8 
I IV : l c T ,y <lr< In, la, when 8 < (V/V R ) 2 < 25 ’ 

[ V : IrCln, la, lT y, when (V/V R ) 2 » 25 

(36) 

which are labelled by the red Roman numbers with the 
boundaries indicated by the red crosses at the upper 
frame of Fig. [3](c). Specifically, when (V/Vr ) 2 < 1, the 
transverse spin diffusion length decreases with the in¬ 
crease of the disorder strength. When 1 < (V/Vr ) 2 < 8 
[(V/Vr ) 2 25], it is obtained from Table IH1 that ft 


Vsinlr/C^la) (Lj, ft \JV3lalr) will 1 the steady-state 
spin polarization being single exponential (oscillation) 
decay. Therefore, the transverse spin diffusion length de¬ 
creases with the increase of the disorder strength. Fur¬ 
thermore, it can be seen in Fig. (3J]c) that the results 
calculated from the analytical formulas Eqs. (IA19D and 
(IA22I) agree with the numerical ones in Regime V. 

For the longitudinal spin diffusion along the y- 
direction, it is shown in Fig.[3|d) that with the increase of 
the scattering strength, for both the strong (450 E r ) and 
weak (10 E r ) Zeeman fields, the longitudinal spin diffu¬ 
sion length is suppressed first and then become insensi¬ 
tive to the scattering. Specifically, when il = 10 E r (l a -C 
In), the longitudinal spin diffusion is robust against the 
scattering except with extremely weak scattering. Here, 
for the strong (weak) Zeeman field H = 450 E r (10 E r ), 
the regimes are divided according to Eq. (IM1) [Eq. (1551) ]. 
For the strong Zeeman field (450£V), when (V/Vr) 2 < 8, 
the longitudinal spin diffusion is suppressed by the scat¬ 
tering. When 8 < (V/V R ) 2 < 50 [(' V/V R f » 50], it is 
obtained from Table HX1 that L v l ft l T l a /(V^ln) (L v l ft l a ) 
with the steady-sate spin polarization being single ex¬ 
ponential decay. Hence, the longitudinal spin diffusion 
length decreases first and then become insensitive to 
the scattering with the increase of the disorder strength. 
For the weak Zeeman field (10.E).), when (V/Vr) 2 < 1, 
our calculation shows that the longitudinal spin diffusion 
length decreases slowly with the increase of the disorder 
strength. When (V/Vr) 2 > 1, one obtains from Table HD 
that L v l ft l a (single exponential decay). Hence, the lon¬ 
gitudinal spin diffusion length is insensitive to the scat¬ 
tering in a wide range. Moreover, it is shown in Fig. 0(d) 
that the results calculated from the analytical formula 
Eq. (IA15D agree with the numerical ones in both the mod¬ 
erate and strong scattering regimes. 

Finally, we address the specific features in the scatter¬ 
ing strength dependence of the transverse and longitu¬ 
dinal spin diffusions along the x- and y-directions. Our 
calculations show that in the weak (l T > In, l a ) scattering 
regime (Regime I), both the transverse and longitudinal 
spin diffusions are suppressed by the scattering. In the 
strong scattering limit (l T <C ln,l a ,l%), the longitudinal 
spin diffusions along the x and y-directions are insen¬ 
sitive to the scattering. Specifically, when l a <C In, the 
longitudinal spin diffusion along the y-direction is robust 
against the scattering except with extremely weak scat¬ 
tering. We emphasize that this robust spin diffusion can¬ 
not be obtained from the over-simplified drift-diffusion 
model, where l s = \JDt s (\l) with D = VpTk/3~V— From 
the drift-diffusion model, in Regime V with r s (k) oc t // 1 , 
one surely obtains that the spin diffusion length is irrel¬ 
evant to the scattering. However, in Regimes I and II, 
i.e., the weak scattering regime defined in the conven¬ 
tional DP spin relaxation ) 38 ! 40 ' 44 r s (k) ft t*,. Hence it 
is obtained that l s oc rfc with the spin diffusion length 
suppressed by the scattering. One further notes that the 
experimental condition can be easily realized to observe 
this robust spin diffusion. On one hand, the spin diffu- 
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sion is set to be along the y-direction with the initial spin 
polarization parallel to the Zeeman field; on the other 
hand, the Zeeman field is tuned to be much weaker than 
the spin-orbit coupled one by the laser field. 


B. Zeeman field dependence 

In this part, we address the Zeeman field dependence 
of the transverse and longitudinal spin diffusions along 
the x- [Fig. 0fa)] and y-directions [Fig. Hfb)] . The SOC 
strength is set to be «o and the spin polarization is chosen 
to be P = 20%. Here, we mainly address the specific 
features in the Zeeman field dependence when the system 
lies in the moderate and strong scattering regimes, which 
can be realized by setting ( V/Vr ) 2 = 60 and ( V/Vr ) 2 = 
10, respectively. 



n/(iooE r ) 

FIG. 4: (Color online) Zeeman field dependence of the trans¬ 
verse and longitudinal spin diffusions along the (a) x- and 
(b) y-directions. The SOC strength a = ao and the spin po¬ 
larization P = 20%. Both the situations with the scattering 
strength (V/Vr) 2 = 60 (squares) and (V/Vr) 2 = 10 (circles) 
calculated numerically are presented. The green (blue) Ro¬ 
man numbers at the upper (lower) frame represent the differ¬ 
ent regimes when (V/Vr) 2 = 10 for the spin diffusion along 
the x-direction (y-direction) with the boundaries indicated 
by the green (blue) crosses. 


We first focus on the case with (V/Vr ) 2 = 60. When 
(V/Vr ) 2 = 60, one observes from Figs. Ufa) and (b) that 
for the spin diffusion along both x- and y-directions, the 
transverse diffusion length decreases with the increase 
of the Zeeman field, as shown by the red solid curve 
with squares. This is because when (V/Vr ) 2 = 60, 
for the spin diffusion along the x-direction (y-direction), 
It < l a ,In, I)/,x (4 < lajnj^y) is satisfied. Hence, the 
system lies in Regime V. Accordingly, for the transverse 
spin diffusion along the x-direction (y-direction), one ob¬ 
tains from Table. [T] (Table. HH that Lj. ft 2l T ln/\/3 

[Lj, ft %/3ZnZ T /(2Z a )]. Therefore, with the increase of 
the Zeeman field, the transverse spin diffusion length de¬ 
creases. For the longitudinal spin diffusion, it is shown 
in Figs. Ufa) and (b) that no matter the spin diffusion is 
along the x-direction or the y-direction, the spin diffu¬ 
sion length is marginally influenced by the Zeeman field 
(blue dashed curve with squares). This anomalous be¬ 
havior can be well understood from the analytical results 
in Regime V. For the longitudinal spin diffusion along 
the x-direction (y-direction), it is obtained that L x l ft l a 
(L v l ft l a ). Accordingly, in Regime V, the longitudi¬ 
nal spin diffusions along both the x- and y-directions 
are marginally influenced by the Zeeman field. Based on 
the drift-diffusion model^ - — it is emphasized that this 
unique feature arises from the insensitivity of diffusion 
coefficient and SRT to the Zeeman field in the strong 
scattering limit. 

We then analyze the case with (V/Vr ) 2 = 10. We first 
divide the regimes for the spin diffusion along the x- and 
y-directions, respectively. For the spin diffusion along 
the x-direction, the regimes for the system are divided 
into 

{ IV : l C T, x <l T < In, l a , when H < 30 E r 

V : l T < In, la,l c T , x , when 30F r < H < 80 E r , 

II : In < l T < l a , when Q > 80 E r 

(37) 

which are labelled by the red Roman numbers with the 
boundaries indicated by the red crosses at the upper 
frame of Fig. |4|a). For the spin diffusion along the y- 
direction, when B < 80 E r and H > 80 E r , the system lies 
in Regimes IV and II, which are represented by the blue 
Roman numbers with the boundaries indicated by the 
blue crosses at the lower frame of Fig.@|b). It is shown in 
Figs.[U[a) and (b) that in Regime II, the transverse (lon¬ 
gitudinal) diffusion length decreases slowly (increases) 
with the increase of the Zeeman field, represented by the 
red solid (blue dashed) curve with circles. This is because 
for the transverse spin diffusion along the x-direction (in¬ 
direction), in Regime II, Lj, ft 2 Z t /i/3 (Lj, ft y/2l T ), with 
the diffusion length marginally influenced by the Zeeman 
field. For the longitudinal spin diffusion, in Regime II, 
L x l ft l T l a /(Vsln) [• L v l ft l T l a /(V3ln)\, leading to the 
enhancement of the spin diffusion by the Zeeman field. 
We emphasize that this enhancement of the spin diffu¬ 
sion arises from the suppression of the longitudinal spin 
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relaxation by the Zeeman field. 

Finally, we emphasize the unique features in the Zee- 
man field dependence of the spin diffusions along the x- 
and y-directions, which can be observed in the experi¬ 
ment. These unique features arise in the longitudinal sit¬ 
uation, i.e., the initial spin polarization is parallel to the 
Zeeman field, for the spin diffusion along both the x- and 
y-directions. On one hand, when the scattering is strong 
with the system in Regime V, the longitudinal spin diffu¬ 
sion is marginally influenced by the Zeeman field; on the 
other hand, when the scattering is relatively weak and 
the Zeeman field is strong with the system in Regime 
II, the spin diffusion length is enhanced by the Zeeman 
field. These unique features can be understood in the 
framework of modified drift-diffusion model addressed in 
Sec.ITTTl 4 - - — It is emphasized that here the drift-diffusion 
model is applicable, which is very different from the lon¬ 
gitudinal spin diffusion along the y-direction in Regimes 
III and IV [the red solid curve with circles in Fig. [3jd)] 
where the modified inhomogeneous broadening picture is 
used. We readdress that with the strong (weak) Zeeman 
and weak (strong) spin-orbit coupled fields, the condition 
to apply the modified drift-diffusion model/modified in¬ 
homogeneous broadening picure is that the initial spin 
polarization is parallel/perpendicular to the larger field. 


C. SOC strength dependence 

In this part, we analyze the SOC strength dependence 
of the transverse and longitudinal spin diffusions along 
the x- [Fig. [5](a)] and y-directions [Fig. 0(b)] . It has 
been well understood from the drift-diffusion model that 
in the strong scattering regime, the spin diffusion length 
is suppressed by the SOC (l s oc l/|ai|)4£r— In our study, 
besides the suppression of the spin diffusion length by the 
SOC, we find two unique features in the SOC strength 
dependence which can not be derived from the over¬ 
simplified drift-diffusion model: the spin diffusion length 
can be either marginally influenced or even enhanced by 
the SOC. These features can be realized when the sys¬ 
tem lies in the moderate and strong scattering regimes. 
Accordingly, in our calculation, the Zeeman field is set 
to be SI = 45 E r - the scattering strengths are set to be 
(V/Vr) 2 = 100 and ( V/Vr ) 2 = 20, respectively. 

We first address the first unique feature, i.e., the 
marginal influence of the SOC on the spin diffusion. In 
Figs. El a) and (b), one observes that when ( V/Vr ) 2 = 
100, no matter the spin diffusion is along the x- or y- 
direction in the transverse configuration (the gray dashed 
curve with squares), when a < 4ao, the spin diffu¬ 
sion length is marginally influenced by the SOC. This 
is because when a < 4ao, the system lies in Regime V 

with Lj, ss \J21 t Iq/V3 (A£ ~ \JV31 t Iq) for the trans¬ 
verse spin diffusion along the x-direction (y-direction). 
It is emphasized that this robustness of the spin dif¬ 
fusion to the SOC cannot be obtained from the drift- 



FIG. 5: (Color online) SOC dependence of the transverse 
and longitudinal spin diffusions along the (a) x- and (b) in¬ 
directions. The Zeeman field O = 45 E r and the spin polar¬ 
ization P = 20%. Both the situations with the scattering 
strength (V/Vr ) 2 = 100 (squares) and (V/Vr ) 2 = 20 (cir¬ 
cles) calculated numerically are presented. The red (blue) 
Roman numbers at the upper (lower) frame represent the dif¬ 
ferent regimes for the spin diffusion along the x-direction (y- 
direction) when (V/Vr ) 2 = 20 with the boundaries indicated 
by the red (blue) crosses. 

diffusion model4£~— As we have addressed in Sec. nm in 
the transverse spin diffusion in Regime V, the inhomo¬ 
geneous broadening has dominant influence on the spin 
diffusion. 

We then analyze the second unique feature where the 
spin diffusion length can be enhanced by the SOC. For 
the spin diffusion along the x-direction, it is shown in 
Fig. El a) by the green (blue) dashed curve with circles 
that when (V/Vr) 2 = 20 with a > 2ao, the transverse 
(longitudinal) spin diffusion is significantly enhanced by 
the SOC; for the spin diffusion along the y-direction, 
when (V/Vr) 2 = 20 with a > 2a 0 the transverse spin 
diffusion length also increases with the increase of the 
SOC (the green dashed curve with circles). This can be 
understood as follows. When a > 2ag, the system lies in 
Regime III. Accordingly, for the transverse (longitudinal) 
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spin diffusion along the x-direction, one obtains L « 
Irln/iVSla) [if, ~ V2l T ln/{V3la)}', for the transverse 
spin diffusion along the y-direction, « y/3l T ln/(2l a ). 
This unique feature is in contrast to the prediction of the 
drift-diffusion modeler— 

In above sections, we have compared the analytical re¬ 
sults [Eqs. (1A91) . (I A 10 1) . (1A15I) . (I A 191) and (1A22[) ] with the 
numerical ones. Now, we address the general condition 
that the analytical results can be applied. It is noted that 
our analytical results are derived in the strong scattering 
regime with l T < l a ,ln and then extended to the moder¬ 
ate scattering regime. The numerical calculations show 
that in the relatively strong and strong scattering regimes 
(l T < Ion l n), the analytical results agree with the numer¬ 
ical ones fairly well; in the moderate scattering regime, 
the condition to use the analytical results is In < l T < la. 
or l a fS It < In-, which are close to the boundary between 
the moderate and relatively strong scattering regimes. 
However, even when the system is away from the bound¬ 
ary between the moderate and strong scattering regimes, 
the dependencies of the spin diffusion on the scattering 
strength, Zeeman field and SOC strength are qualita¬ 
tively correct in the moderate scattering regimes. 


V. CONCLUSION AND DISCUSSION 

In conclusion, we have investigated the steady-state 
spin diffusion for the 3D ultracold spin-orbit coupled 
40 K gas by the KSBE approach^ 4 first analytically and 
then numerically. The spin diffusions along the x- and 
y-directions for the transverse (P||z) and longitudinal 
(P| |x) configurations are studied. It is first shown analyt¬ 
ically that the behaviors of the steady-state spin diffusion 
in the four configurations (x-T, x-L, y-T and y-L) are 
determined by three characteristic lengths: the mean free 
path l T , the Zeeman oscillation length In, and the SOC 
oscillation length l a . We have derived the spin diffusion 
lengths for the spin diffusions in the four configurations in 
the strong scattering regime, which are then extended to 
the weak scattering one. We further find that in different 
limits, the complex analytical reuslts can be reduced to 
different extremely simple forms, and correspodingly, the 
system can be divided into different regimes. Specifically, 
it is revealed that by tuning the scattering strength, the 
system can be divided into five regimes: I, weak scat¬ 
tering regime (l T > In, la)', II, Zeeman field-dominated 
moderate scattering regime (In <C l T <C l a )', III, SOC- 
dominated moderate scattering regime (l a -C l T <C In)', 

IV, relatively strong scattering regime ( l % < Z T < ln,l a )', 

V, strong scattering regime (l T <C In, la, It)- I n different 
regimes, the corresponding behaviors of the spatial evo¬ 
lution of the spin polarization in the steady state are 
very rich, showing different dependencies on the scatter¬ 
ing strength, Zeeman field and SOC strength. These de¬ 
pendencies are summarized in Table U (Table [[]} for the 
spin diffusion along the x-direction (y-direction). Then 
the scattering strength, Zeeman field and SOC strength 


dependencies of the spin diffusions are numerically cal¬ 
culated and compared with the analytical ones. It is 
shown that the analytical results agree with the numeri¬ 
cal ones fairly well in the relatively strong/strong scatter¬ 
ing regime and the region close to the boundary between 
the moderate and relatively strong scattering regimes. 
However, it is found that even when the system is away 
from the strong scattering regime, the analytical results 
are still qualitatively correct in the moderate scattering 
regimes. 

The rich behaviors of the spin diffusions in differ¬ 
ent regimes are hard to be understood in the frame¬ 
work of the previous simple drift-diffusion model^— 0 or 
the direct inhomogeneous broadening [Eqs. © and ©] 
picture 44 *^ 54 in the literature. In this work, we extend 
our previous inhomogeneous broadenings [Eqs. © and 
©] to the effective ones in Eqs. (1201) and (12Tj) [Eqs. (1311) 
and (p© [ for the spin diffusion along the x-direction (in¬ 
direction). In the limit situation, we suggest reasonable 
pictures referred to as modified drift-diffusion model and 
modified inhomogeneous broadening picture to facilitate 
the understanding of the simple analytical results in Ta¬ 
bles [I] and [© It is shown that the behaviors of the spin dif¬ 
fusions can be analyzed in the situation either with strong 
Zeeman and weak spin-orbit coupled fields (Regimes II 
and V) or weak Zeeman and strong spin-orbit coupled 
fields (Regimes III and IV). When the spin polarization 
is parallel (perpendicular) to the larger field between the 
Zeeman and spin-orbit coupled fields, the spin polariza¬ 
tion cannot (can) rotate around the effective inhomoge¬ 
neous broadening fields efficiently, and hence the modified 
drift-diffusion model (modified inhomogeneous broaden¬ 
ing picture) can be used. In the modified drift-diffusion 
model, in the strong scattering regime, t s ( k) remains 
the SRT in the conventional DP mechanism j 29 i 30 i 32 i 44 
whereas in the moderate scattering regime, r s (k) is re¬ 
placed by the helix spin-flip rates determined in the he¬ 
lix space.- 1,39 In the modified inhomogeneous broaden¬ 
ing picture, the behavior of the spin diffusion is deter¬ 
mined by the effective inhomogeneous broadenings to¬ 
gether with the spin-conserving scattering. Based on the 
modified drift-diffusion model and modified inhomoge¬ 
neous broadening picture, apart from Regime IV, all the 
features in different regimes can be obtained. Below, we 
address several anomalous features of the spin diffusion, 
which are in contrast to both the simple drift-diffusion 
model and the direct inhomogeneous broadening picture. 

In the scattering strength dependence, it is found that 
when l a <C In, the longitudinal spin diffusion along the 
y-direction is robust against the scattering even when the 
system is away from the strong scattering regime, which 
is in contrast to the simple drift-diffusion model. In that 
model, in the weak scattering regime, with l s oc k^/rf/m, 
the spin diffusion length is suppressed by the scatter¬ 
ing. In the Zeeman field dependence, when the system is 
in Regime II, the longitudinal spin diffusion is enhanced 
by the Zeeman held. This is in contrast to the pre¬ 
diction from the previous drift-diffusion model and the 












16 


direct inhomogeneous broadening picture. In the sim¬ 
ple drift-diffusion model, in the strong (weak) scatter¬ 
ing regime, with l s oc l/(m|a|) ( l s oc k^/rf/m), it is 
obtained that the diffusion length is irrelevant to the 
Zeeman field. In the direct inhomogeneous broadening 
picture, the spin diffusion is suppressed due to the en¬ 
hancement of the inhomogenous broadening by the Zee- 
man field. Finally, in the SOC strength dependence, we 
find that the spin diffusion length can also be enhanced 
by the SOC in Regime III. This also goes beyond the 
prediction from the simple drift-diffusion model and the 
direct inhomogeneous broadening picture. In the sim¬ 
ple drift-diffusion model, in the strong (weak) scattering 
regime, with l s oc l/(m|a|) (l s oc k^/rf/m), the spin dif¬ 
fusion length is suppressed (uninfluenced) by the SOC. 
In the direct inhomogeneous broadenings picture, the 
spin diffusion is uninfluenced (suppressed) for the spin 
diffusion along the x-direction (y-direction). All these 
anomalous behaviors have been well understood from our 
modified drift-diffusion model and/or modified inhomo¬ 
geneous broadening picture. 

We emphasize that for the longitudinal spin diffusion 
along the y-direction, the robustness against the scatter¬ 
ing strength exists in a wide range including both the 
strong and weak scattering regimes. It is noted that 
in the strong scattering regime, this feature has been 
predicted in the simple drift-diffusion model^ - — and is 
also revealed in graphene by the KSBE approach— In 
this work, we further extend it into the weak scattering 
regime. Moreover, it is found that in a wide range of the 
scattering, the corresponding diffusion length is only de¬ 
termined by the SOC strength, and hence also irrelevant 
to the atom density and temperature. One expects sim¬ 
ilar behavior in symmetric (110) quantum wells under a 
weak in-plane magnetic field with similar SOC i 38 i 39 i 64 — 
Moreover, the enhancement of the longitudinal spin dif¬ 
fusion by the Zeeman field has not yet been reported in 
the literature, which is also expected to be observed in 
symmetric (110) quantum wells when the Zeeman energy 
larger than the spin-orbit coupled one. 

Finally, we further compare the pictures of the spin dif¬ 
fusion provided in this work to understand the calculated 
results and those in the literature. In the previous works, 
the spin diffusion in the strong scattering regime have 
been extensively studied in the system with SOC, includ¬ 
ing semiconductors ) 44 ! 69- — graphen o 50 ’ 76- — and recently 
monolayer M 0 S 2 — , Drift-diffusion model^ - — and/or the 
inhomogeneous broadening picture 44 ^ - — were used to 
understand the behaviors of the spin diffusion. In this 
work, the analytical results in the strong scattering 
regime are extended to the weak one and confirmed by 
the full numerical calculation. For the strong scatter¬ 
ing regime, we further divide it into the relatively strong 
and strong scattering regimes; whereas for the weak scat¬ 
tering regime, we divide it into the moderate and weak 
scattering regimes. We find all the anomalous behaviors 
revealed in this work appear in the moderate and rela¬ 
tively strong scattering regimes. Furthermore, our modi¬ 


fied drift-diffusion model and/or modified inhomogeneous 
broadening picture are used to understand the behaviors 
of the spin diffusion in all these regimes. It is found that 
in the moderate and strong scattering regimes, these pic¬ 
tures work well. However, in the relatively strong scat¬ 
tering regime, which lies in the crossover of the moder¬ 
ate and strong scattering regimes, these pictures fail to 
explain the behaviors of the spin diffusion along the y- 
direction. This is because for the drift-diffusion model, 
in the relatively strong scattering regime, it is too rough 
to consider the anisotropy between the diffusions along 
different directions— Nevertheless, when the scattering 
is strong enough, this anisotropy in the spin diffusion 
behavior is vanished. Whereas for the inhomogeneous 
broadening picture, in this regime there exists strong 
competition between the effective inhomogeneous broad¬ 
ening and scattering, which makes the behavior of the 
spin diffusion complicated ! 50 ’ 53 ' 54 
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Appendix A: Analytical Analysis 


We analytically derive the transverse and longitudinal 
spin diffusion lengths for the spin diffusions along the st¬ 
and y-directions based on the KSBEs [Eq. (TH?1) ]. 

Generally, the density matrix depends on both the 
zenith (between k and 0-axis) and azimuth (between k 
and (/-axis in the y-z plane) angles 0 k and 0k in 3D. 
However, with the specific form of the SOC [Eq. (JT])] and 
isotropic scattering terms [Eq. Q], we can define the 
quantity 




Pk — 2^ I ^^kPki 


(Al) 


which is averaged over the azimuth angle 0k, to describe 
the kinetics of the density matrix—. Accordingly, the KS¬ 
BEs in the steady state become 


h d Pk(r) 


flcr x /2,p k (r) 


m dt; 

+ Wkk ' [£k( r ) _ Pk'( r )] = 0, 


ak x cr z /2,p k (r) 

(A2) 


in which pk(r) only depends on 0^. 
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1. Spin diffusion along the x-direction 


2. Spin diffusion along the y-direction 


For the spin diffusion along the x-direction, pk is ex¬ 
panded by the Legendre function, which is written as 


Pk = 'y'p[C?Pi{cos0 k ), 


l 


(A3) 


with Cf = \/(2l + 1)/(47t). Accordingly, the dynamical 
equation for p\ is written as Eq. ffTTTli . With the spin 
vector defined by S l k = Tr [p l k a] , the equations for the 
spin vectors can be obtained. By further keeping the 
zeroth and first orders (l = 0,1), the equation for the 
vector S fc = (SjJ iX , S^, S^, S^, S^, S£J T is written 
as 


J k,xi ^k,y’> ^k : zi VJ k,y'> ^k,z) 

d x S k + U x S k = 0, 


with 


U x = 


/ 

0 1/4 0 

Vs/ir 

0 

0 


-1/4 

0 0 

0 

VS/l T 

1/4 


0 

0 0 

0 

-1/4 

VS/l T 


0 

0 0 

0 

i/4 

0 


0 

0 1/4 

-1/4 

0 

0 

V 

0 -1/4 o 

0 

0 

0 
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(A5) 

From Eq. (IA5I) . the spin diffusion and oscillation lengths 
can be found from the eigenvalues of U x denoted by X x , 
which satisfy 


X x H- a\ x bX x -\- c — 0, 


(A6) 


with a = 2/Iq + 2/Z^,, 6 = 3/(44) 2 + 1 /I q + 2 /(44) 2 + 
1 /l% and c = —3/(444) 2 - In Eq. (IA6I) . the real 
and imaginary parts of 1/X x correspond to the diffusion 
length and oscillation length, respectively. 

With A = (q/2) 2 + (p/3) 3 where q = 2a 3 /27 — ab/S + c 
and p = b — a 2 / 3, it is demonstrated that in the mod¬ 
erate and strong scattering regimes, A is always larger 
than zero. Therefore, there are one real root (A re ) and 
two complex conjugate roots (Ai mj ±) for Af. in Eq. 
which are written as 


A re = \]-q/2-VA + ^J-q/2 + VA, (A7) 

A im ,± = -(l/2)^-q/2-VA + \j—q/2 + Va) 

± (VSi/2) ~q/2 -VK- \J ~q/2 + Va) . (A8) 

Accordingly, the spin diffusion length for the single expo¬ 
nential decay, the spin diffusion length for the oscillation 
decay and the spin oscillation length are given by 

L x s =1/v^_ (A9) 

L x 0 = v^/yAr+^/A^^~7^A~Tp/3, (A10) 

^ = 2v / 3LS/|Aim, + -A im ,_|. (All) 


For the spin diffusion along the y-direction, pk is ex¬ 
panded by the Fourier function, which is denoted by 


P k 


= ^2Pk^p{ilO k ). 


(A12) 


The corresponding dynamical equation for p k is writ¬ 
ten as Eq. (l22l) . By keeping the zeroth and first orders 
(l = 0,1), the equation for the vector S*, = (S° x , S/ x — 
Sfc, 1 *, Sg >w , S° >z , Sl iV - S^, S\ <z - S^) T is written as 


dyS k + UyS k = 0, 


(A4) where 


/ 


Uy — 


ivmi+.uL 


\ 


0 

0 

0 

0 

il T 

~T 

l OL 

2 i 


0 

0 

0 

0 - 
2 i 


(A13) 

t \ 


l 

It V 
i i 


Vsk 


Vsk 


Vsin 

o 


Vs 4 

0 


It 


0 


/ 

(A14) 

In above equation, it is noted that the up-left 2x2 and 
down-right 4x4 blocks C/“ _1 and Uy ~ r , which describe the 
spin diffusion for and S y /S-, are decoupled to each 
other. 

Accordingly, from the eigenvalues of I/”' 1 , the spin dif¬ 
fusion length for S x is found to be 

L v L = l a yJl 2 T /(3l 2 n) + l. (A15) 

From Uy~ r , it is found that the four eigenvalues X y satisfy 
[4/(3 Zq) — 1/Zq] X 2 +4/(31q1 2 )+4/(91q) = 0. (A16) 


From Eq. (|A16[) . when 1/Z^ — 8/(SIq)(1/1^ + 2/1 2 ) > 0, 
there are four real roots, among which the two positive 
ones read 


\K\ = X 


1 

2Zl 


3 lh 


zb , 




o /2 /2 




(A 17) 


Accordingly, the steady-state spin polarization S y or S 2 
is limited by the bi-exponential decay, with the diffusion 
length being 


l v ^ + = i/|a; 

L y T~ = 1/| A. 


(A18) 

(A19) 
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respectively. Otherwise, when 1 /l* — 8/(3Iq)(1// 2 + 
2 /l\) < 0, the four roots for X y are complex. Specifi¬ 
cally, the real part of the roots for Eq. (1A16I) is identical, 
which is written as 

A r ; = \J 1/(4Z2) - 17(3*2) + ^1/(9 + 1/(2,1111). 

(A20) 

The complex part of the roots (absolute value) for 
Eq. (IA16I) is 

|A‘ m | = + (A21) 


Therefore, the steady-state spin polarization S y or S, is 
determined by the oscillation decay with the decay length 
and oscillation length being 


L v t = 1/A™, (A22) 

l y T = 2A™/|A”|, (A23) 


respectively. 


* Author to whom correspondence should be addressed; 
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1 J. Y. Vaishnav and C. W. Clark, Phys. Rev. Lett. 100, 
153002 (2008). 

2 J. Y. Zhang, S. C. Ji, Z. Chen, L. Zhang, Z. D. Du, B. 
Yan, G. S. Pan, B. Zhao, Y. J. Deng, H. Zhai, S. Chen, 
and J. W. Pan, Phys. Rev. Lett. 109, 115301 (2012). 

3 C. Qu, C. Hamner, M. Gong, C. Zhang, and P. Engels, 
Phys. Rev. A 88, 021604(R) (2013). 

4 L. J. LeBlanc, M. C. Beeler, K. J. Garcia, A. R. Perry, S. 
Sugawa, R. A. Williams, and I. B. Spielman, New J. Phys. 
15, 073011 (2013). 

5 M. C. Beeler, R. A. Williams, K. J. Garcia, L. J. LeBlanc, 
A. R. Perry, and I. B. Spielman, Nature 498, 201 (2013). 

6 P. L. Pedersen, M. Gajdacz, F. Deuretzbacher, L. Santos, 
C. Klempt, J. F. Sherson, A. J. Hilliard, and J. J. Arlt, 
Phys. Rev. A 89, 051603(R) (2014). 

7 Y. J. Lin, K. J. Garcia, and I. B. Spielman, Nature (Lon¬ 
don) 471, 83 (2011). 

8 K. J. Garcia, L. J. LeBlanc, R. A. Williams, M. C. Beeler, 
C. Qu, M. Gong, C. Zhang, and I. B. Spielman, Phys. Rev. 
Lett. 114, 125301 (2015). 

9 A. Sommer, M. Ku, and M. W. Zwierlein, New J. Phys. 
13, 055009 (2011). 

10 H. Heiselberg, Phys. Rev. Lett. 108, 245303 (2012). 

11 O. Goulko, F. Chevy, and C. Lobo, Phys. Rev. Lett. Ill, 
190402 (2013). 

12 C. Cao, E. Elliott, J. Joseph, H. Wu, J. Petricka, T. 
Schafer, and J. E. Thomas, Science 331, 58 (2011). 

13 A. Sommer, M. Ku, G. Roati, and M. W. Zwierlein, Nature 
472, 201 (2011). 

14 G. M. Bruun, New J. Phys. 13, 035005 (2011). 

15 T. Enss and R. Haussmann, Phys. Rev. Lett. 109, 195303 

( 2012 ). 

16 T. Enss, Phys. Rev. A 88, 033630 (2013). 

11 A. J. Leggett and M. J. Rice, Phys. Rev. Lett. 20, 586 
(1968); A. J. Leggett, J. Phys. C 3, 448 (1970). 

18 M. Koschorreck, D. Pertot, E. Vogt, and M. Kohl, Nat. 
Phys. 9, 405 (2013). 

19 A. B. Bardon, S. Beattie, C. Luciuk, W. Cairncross, D. 
Fine, N. S. Cheng, G. J. A. Edge, E. Taylor, S. Zhang, S. 
Trotzky, and J. H. Thywissen, Science 344, 722 (2014). 

20 S. Trotzky, S. Beattie, C. Luciuk, S. Smale, A. B. Bardon, 
T. Enss, E. Taylor, S. Zhang, and J. H. Thywissen, Phys. 
Rev. Lett. 114, 015301 (2015). 

21 T. Enss, Phys. Rev. A 91, 023614 (2015). 


22 X. Du, L. Luo, B. Clancy, and J. E. Thomas, Phys. Rev. 
Lett. 101 , 150401 (2008). 

23 F. Piechon, J. N. Fuchs, and F. Laloe, Phys. Rev. Lett. 
102, 215301 (2009). 

24 S. S. Natu and E. J. Mueller, Phys. Rev. A 79, 051601(R) 
(2009). 

25 X. Du, Y. Zhang, J. Petricka, and J. E. Thomas, Phys. 
Rev. Lett. 103, 010401 (2009). 

26 P. Wang, Z. Q. Yu, Z. Fu, J. Miao, L. Huang, S. Chai, H. 
Zhai, and J. Zhang, Phys. Rev. Lett. 109, 095301 (2012). 

2| L. W. Cheuk, A. T. Sommer, Z. Hadzibabic, T. Yefsah, 
W. S. Bakr, and M. W. Zwierlein, Phys. Rev. Lett. 109, 
095302 (2012). 

28 U. Ebling, J. S. Krauser, N. Flaschner, K. Sengstock, C. 
Becker, M. Lewenstein, and A. Eckardt, Phys. Rev. X 4, 
021011 (2014). 

29 I. V. Tokatly and E. Ya. Sherman, Phys. Rev. A 87, 
041602(R) (2013). 

30 S. S. Natu and S. Das Sarma, Phys. Rev. A 88, 033613 
(2013). 

31 T. Yu and M. W. Wu, Phys. Rev. A 88, 043634 (2013). 

32 J. Radic, S. S. Natu, and V. Galitski, Phys. Rev. Lett. 
112, 095302 (2014). 

33 I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

34 S. S. Kondov, W. R. McGehee, J. J. Zirbel, and B. De¬ 
Marco, Science 334, 66 (2011). 

35 F. Jendrzejewski, A. Bernard, K. Miiller, P. Cheinet, V. 
Josse, M. Piraud, L. Pezze, L. S. Palencia, A. Aspect, and 
P. Bouyer, Nat. Phys. 8, 398 (2012). 

36 D. Clement, A. F. Varon, J. A. Retter, L. S. Palencia, A. 
Aspect, and P. Bouyer, New J. Phys. 8, 165 (2006). 

37 M. Piraud, L. Pezze, and L. S. Palencia, New J. Phys. 15, 
075007 (2013). 

38 Y. Zhou, T. Yu, and M. W. Wu, Phys. Rev. B 87, 245304 
(2013). 

39 T. Yu and M. W. Wu, Phys. Rev. B 89, 045303 (2014). 

40 M. I. D’yakonov and V. I. Perel’, Zh. Eksp. Teor. Fiz. 60, 
1954 (1971) [Sov. Phys. JETP 33, 1053 (1971)]. 

41 Y. Yafet, Phys. Rev. 85, 478 (1952). 

42 R. J. Elliott, Phys. Rev. 96, 266 (1954). 

43 J. P. Brantut, J. Meineke, D. Stadler, S. Krinner, and T. 
Esslinger, Science 31, 337 (2012). 

44 M. W. Wu, J. H. Jiang, and M. Q. Weng, Phys. Rep. 493, 
61 (2010). 

45 M. Ziese and M. J. Thornton, eds., Spin Electronics 










19 


(Springer, Berlin, 2001). 

46 Z. G. Yu and M. E. Flatte , Phys. Rev. B 66, 201202(R) 

( 2002 ). 

4 ' I. Zutic , J. Fabian, and S. Das Sarma, Phys. Rev. Lett. 
88, 066603 (2002). 

48 J. Fabian, I. Zutic , and S. Das Sarma, Phys. Rev. B 66, 
165301 (2002). 

49 J. Fabian, A. M. Abiague, C. Ertler, P. Stano, and I. Zutic, 
Acta Phys. Slov. 57, 565 (2007). 

50 P. Zhang and M. W. Wu, Phys. Rev. B 84, 045304 (2011). 

51 M. W. Wu and C. Z. Ning, Eur. Phys. J. B. 18, 373 (2000); 
M. W. Wu, J. Phys. Soc. Jpn. 70, 2195 (2001). 

52 M. Q. Weng and M. W. Wu, Phys. Rev. B 66, 235109 
(2002); J. Appl. Phys. 93, 410 (2003). 

53 J. L. Cheng and M. W. Wu, J. Appl. Phys. 101, 073702 
(2007). 

54 J. L. Cheng, M. W. Wu, and I. C. da Cunha Lima, Phys. 
Rev. B 75, 205328 (2007). 

55 C. A. Regal, C. Ticknor, J. L. Bohn, and D. S. Jin, Phys. 
Rev. Lett. 90, 053201 (2003). 

56 M. Egorov, B. Opanchuk, P. Drummond, B. V. Hall, P. 
Hannaford, and A. I. Sidorov, Phys. Rev. A 87, 053614 
(2013). 

57 T. Ozawa, L. P. Pitaevskii, and S. Stringari, Phys. Rev. A 
87, 063610 (2013). 

58 R. A. Williams, M. C. Beeler, L. J. LeBlanc, K. J. Garcia, 
and I. B. Spielman, Phys. Rev. Lett. Ill, 095301 (2013). 

59 M. W. Wu and H. Metiu, Phys. Rev. B 61, 2945 (2000). 

60 H. Haug and A. P. Jauho, Quantum Kinetics in Transport 
and Optics of Semiconductors (Springer, Berlin, 1996). 

61 M. Q. Weng and M. W. Wu, Phys. Rev. B 68, 075312 
(2003). 

62 J. H. Jiang and M. W. Wu, Phys. Rev. B 79, 125206 
(2009). 

63 A. F. Ioffe and A. R. Regel, Prog. Semicond. 4, 237 (1960). 

64 Y. Ohno, R. Terauchi, T. Adachi, F. Matsukura, and H. 
Ohno, Phys. Rev. Lett. 83, 4196 (1999); Physica E 6, 817 
(2000); T. Adachi, Y. Ohno, F. Matsukura, and H. Ohno, 
Pliysica E 10, 36 (2001). 

65 M. W. Wu and M. Kuwata-Gonokami, Solid State Com- 


mun. 121, 509 (2002). 

66 S. Dohrmann, D. Hagele, J. Rudolph, M. Bichler, D. 
Schuh, and M. Oestreich, Phys. Rev. Lett. 93, 147405 
(2004). 

67 G. M. Muller, M. Romer, D. Schuh, W. Wegscheider, J. 
Hiibner, and M. Oestreich, Phys. Rev. Lett. 101, 206601 
(2008). 

68 I. V. Tokatly and E. Y. Sherman, Phys. Rev. B 82, 
161305(R) (2010). 

69 T. Ando, A. B. Fowler, and F. Stern, Rev. Mod. Phys. 54, 
437 (1982). 

70 F. Meier and B. P. Zakharchenya, Optical Orientation 
(North-Holland, Amsterdam, 1984). 

71 Semiconductor Spintronics and Quantum Computation, 
edited by D. D. Awschalom, D. Loss, and N. Samarth 
(Springer, Berlin, 2002). 

12 I. Zutic, J. Fabian, and S. D. Sarma, Rev. Mod. Phys. 76, 
323 (2004). 

' 3 J. Fabian, A. Matos-Abiague, C. Ertler, P. Stano, and I. 

Zutic, Acta Phys. Slov. 57, 565 (2007). 

' 4 Spin Physics in Semiconductors, edited by M. I. D’yakonov 
(Springer, Berlin, 2008). 

75 T. Korn, Phys. Rep. 494, 415 (2010). 

' 6 N. Tombros, C. Jozsa, M. Popinciuc, H. T. Jonkman, and 
B. J. van Wees, Nature 448, 571 (2007). 

" C. Ertler, S. Konschuh, M. Gmitra, and J. Fabian, Phys. 
Rev. B 80, 041405(R) (2009). 

78 K. Pi, W. Han, K. M. McCreary, A. G. Swartz, Y. Li, and 
R. K. Kawakami, Phys. Rev. Lett. 104, 187201 (2010). 

' 9 M. H. D. Guimaraes, A. Veligura, P. J. Zomer, T. Maassen, 
I. J. Vera-Marun, N. Tombros, and B. J. van Wees, Nano 
Lett. 2012, 3512 (2012). 

80 W. Han, K. M. McCreary, K. Pi, W. H. Wang, Y. Li, H. 
Wen, J. R. Chen, and R. K. Kawakami, J. Magn. Magn. 
Mater. 324, 369 (2012). 

81 A. Dankert, M. V. Kamalakar, J. Bergsten, and S. P. Dash, 
Appl. Phys. Lett. 104, 192403 (2014). 

82 L. Wang and M. W. Wu, Phys. Rev. B 89, 205401 (2014). 



